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Abstract 

If the X(3872) is a loosely-bound molecule of the charm mesons D°D*° and D*°D°, it can decay 
through the decay of a constituent in a hadronic channel with a nearby threshold, such as J/ipu 
or J/ipp. The differential decay rates of the X into J/ip7r + TT~ , J/ipir + 7t~ir°, J/ipTr 0, y, and J/ipj 
are calculated in terms of XJ/ipp and XJ/iJjuj coupling constants using an effective lagrangian 
that reproduces the decay rates of the u and the p. The dependence of the coupling constants on 
the binding energy and the total width of the X is determined by a factorization formula. Results 
from a model by Swanson are used to predict the partial width of X into J /if) Tr + ir~ir° as a function 
of the binding energy and the total width of the X. 

PACS numbers: 12.38.-t, 12.39.St, 13.20. Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a narrow resonance near 3872 MeV discovered by the Belle collaboration 
in electron-positron collisions through the £?-meson decay B ± — > XK ± followed by the decay 
X — > J/ip7T + 7T~ [1]. Its existence has been confirmed by the CDF and D0 collaborations 
through its inclusive production in proton-antiproton collisions [2, 3] and by the Babar 
collaboration through the discovery mode B ± — > XK ± [4]. The combined measurement of 
the mass of the X is [5] 

m x — 3871.9 ± 0.5 MeV, (1) 

which is within 1 MeV of the threshold for the charm mesons D° and D*°. The presence of the 
J/ijj among the decay products of the X(3872) motivates its interpretation as a charmonium 
state with constituents cc [5-9] . Two possibilities motivated by the proximity of the mass in 
Eq. (1) to the D°D*° threshold are a hadronic molecule with constituents DD* [10-20] and a 
"cusp" at the D°D*° threshold associated with strong coupling to D°D*° or D*°D° [21, 22]. 
Other proposed interpretations include a tetraquark with constituents ccqq [23], a "hybrid 
charmonium" state with constituents ccg [24, 25], a glueball with constituents ggg [26], and 
a diquark-antidiquark bound state with constituents cu + cu [27]. The interpretation as a 
DD* molecule is particularly predictive because the small binding energy implies that the 
molecule has universal properties that are completely determined by the binding energy 
[11, 13, 15, 18]. The small binding energy can be further exploited through factorization 
formulas for production and decay rates of the X [19]. 

Measurements of the decays of the X can be used to determine its quantum numbers and 
narrow down the possibilities [5, 28-32]. The upper bound on the decay width of the X is 

[1] 

T x < 2.3 MeV (90% C.L.), (2) 

which is much narrower than other charmonium states above the DD threshold. The product 
of the branching fractions associated with the discovery channel is [1, 4, 5] 

Br[£+ -> XK + ] Br[X -> J/^tt+tt-] = (1.3 ± 0.3) x 10~ 5 . (3) 

The invariant mass distribution of the two pions from the decay X — > J/ip-K + -K~ seems 
to peak near the upper endpoint, which suggests that the pions come from a virtual p 
resonance. Recently the Belle collaboration has observed the X(3872) in the decay mode 
X — > J/ip 7r + 7r~7T° [32] with the branching ratio 

Br[X -> JMtt+tt-tt ] 

— [ — - J = 1.0 ± 0.4 stat ± 0.3 syst . 4 

Br[X — > J/ip-K+n-} 

The invariant mass distribution of the three pions indicates that they come predominantly 
from a virtual uj resonance [33] . If the decays into J/ ifj 7r + 7r~ and J/ip 7r + 7r~7r° are interpreted 
as J/ip p* and J/ipu*, the approximate equality of the branching fractions in Eq. (4) implies 
a large violation of isospin symmetry. The Belle collaboration also reported evidence for the 
decay X — > [32] with the branching ratio 

Mx^JHn] , 014±0 . 05 . (5) 



Bl[X -> J/iP 71 + 71- 
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The observation of the decay into J/ipl establishes the charge conjugation of the X to be +. 
By analyzing angular distributions in the decay of X into J/ip the Belle collaboration 

has ruled out all J p+ assignments for X with J < 2 other than 1 ++ and 2 ++ [32]. Upper 
limits have been placed on the branching fractions for other decay modes of the X, including 
D°D°, D+D-, D°D 7i° [34], X dl, Xc2l, J/V^ ^ [35], and J/lpr] [36]. Upper limits have 
also been placed on the partial widths for the decay of X into e + e~ [37, 38] and into 77 [37]. 

The possibility that charm mesons might form molecular states was considered shortly 
after the discovery of charm [39-42]. In 1993, Tornqvist made a quantitative study of the 
possibility of molecular states of charm mesons using a one-pion-exchange potential model 
[43]. He found that the isospin-0 combinations of DD* and D*D could form weakly-bound 
states in the S-wave 1 ++ channel and in the P-wave h channel. After the discovery of the 
X(3872), Tornqvist pointed out that because the binding energy is small compared to the 
splitting between the D + D*~ and D°D*° thresholds, there will be large violations of isospin 
symmetry [10] . Swanson considered a potential model that includes both one-pion-exchange 
and quark exchange and found that the C = + superposition of D°D*° and D*°D° could 
form a weakly-bound state in the S-wave 1 ++ channel [14]. Swanson's model included not 
only the charmed mesons D°D*°, D*°D°, D + D*~, and D* + D~, but also two other pairs 
of hadrons with nearby thresholds: J/ip p° and J/ipu. His prediction that the branching 
fraction for the decay of X into J/ip 7r + 7r~7r° should be comparable to that for decay into 
J/ipir + 7T~ was verified by the Belle collaboration [32, 33]. 

In this paper, we analyze decays of the X(3872) into J/ip and light hadrons under the 
assumption that X is a loosely-bound DD* molecule and that these decays proceed through 
transitions oi X to J/ip p and J/ipu. In Section II, we summarize some of the universal 
results for a system with large scattering length and we give the current constraints on the 
real and imaginary parts of the large scattering length for the DD* system. In Section III, we 
discuss other hadronic states with thresholds near the D°D*° threshold and we summarize 
results from Swanson's model for the X(3872). In Section IV, we calculate the differential 
distributions for decays of X into J/ip tt + tt~, J/ipii + 7i^7T , and J/ip7r 0r y using an effective 
lagrangian that reproduces decays of the vector mesons. We also use vector meson dominance 
to calculate the partial width for the decay into J/ip 7. The normalizations of the decay rates 
are determined by unknown coupling constants for the transitions of X to J/ip p and J/ipuj. 
The dependence of the coupling constants on the binding energy and the total width of the 
X is deduced using a factorization formula. In Section V, we use a two-channel scattering 
model to illustrate how the coupling constants can be related to probabilities for the J/ip p 
and J/ipcu components of the X. We use the probability for J/ipuj in Swanson's model 
to give a quantitative prediction for the partial decay rate of the X into J/ip 7r + 7r~7r° as a 
function of the binding energy and the total width of the X. A summary of our results is 
given in Section VI. An updated determination of the parameters in the effective lagrangian 
for light pseudoscalar and vector mesons is given in an Appendix. 

II. UNIVERSALITY AND THE DD* SYSTEM 

The mass of the X is extremely close to the D°D*° threshold: m^o +m£)»o = 3871.3 ±1.0 
MeV. From the mass measurement in Eq. (1), the difference is 

m x - (m D o +m D *o) = +0.6 ± 1.1 MeV. (6) 
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Most of the uncertainty comes from the experimental uncertainty in 2m £,0, because the 
mass difference m^.o — m^o has a much smaller uncertainty. If the X were a D D*°/D*°D 
molecule, the energy difference in Eq. (6) would have to be negative, corresponding to a 
positive binding energy defined by 

E x = (m D o + m D *o) -m x - (7) 

The measurement of the binding energy in Eq. (6) is compatible with a small negative value 
corresponding to a DD* molecule. However, the central value of the energy difference in 
Eq. (6) is positive, corresponding to a resonance in D°D*° and D*°D° scattering rather than 
a bound state. Bugg has referred to this possibility as a "cusp state" [21, 22], because the 
line shape of the X in some of its decay modes has a cusp at the D*°D° threshold. 

The energy difference in Eq. (6) is tiny compared to the natural energy scale for binding 
by the pion exchange interaction: m^/2/i w 10 MeV, where /i is the reduced mass of D° and 
D*°: 

mr,omn*o , s 

H = — - — = 966.5 ± 0.3 MeV. 8 

m D o + m D ,a 

Whether the energy difference is positive or negative, its unnaturally small value implies 
that if the X couples to D°D*° and D*°D°, the S-wave scattering lengths for those channels 
must be large compared to the natural length scale l/m n associated with the pion exchange 
interaction. Since the experimental evidence favors the charge conjugation quantum number 
C = +, we assume that there is a large scattering length a in the C = + channel and that 
the scattering length in the C = — channel is negligible in comparison. In this case, the 
scattering lengths for elastic D°D*° scattering and elastic D*°D° scattering are both a/2. 

Nonrelativistic few-body systems with short-range interactions and a large scattering 
length have universal properties that depend on the scattering length but are otherwise 
insensitive to details at distances small compared to \a\ [44]. We consider the scattering 
length to be large if it is much larger than the natural momentum scale associated with low- 
energy scattering. The universal results are encoded in the truncated connected transition 
amplitude, which is a function of the total energy E of the two particles in the rest frame 
of the pair: 

A(E) = - , (9) 

where \x is the reduced mass of the two particles. If a is real and positive, this amplitude has 
a pole on the physical sheet at E — — l/(2/ia 2 ), indicating the existence of a weakly-bound 
state with the universal binding energy 

E * = ^ < 10 > 
The universal momentum-space wavefunction of this bound state is 

TO = 2 , , / 2 - I 11 ) 

The universal amplitude for transitions from the bound state to a scattering state consisting 
of two particles with small relative momentum is determined by the residue of the pole in 
A(E): 



-4x = — a~ 1/2 . (12) 
/i 
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If there is an inelastic scattering channel, the large scattering length a has a negative 
imaginary part. It is convenient to express the complex scattering length in the form 

-=7re + ^7im, (13) 

a 

where 7 rc and 7; m are real and 7i m > 0. The universal expression for the binding energy in 
Eq. (10) has an imaginary part iTx/2, where 

T X = 27rc7im//i. (14) 

If 7im < 7rc, Tx is the full width at half maximum of a resonance in the inelastic channel 
[19]. It therefore can be interpreted as the rate for the decay of the bound state into the 
inelastic channel. The peak of the resonance is below the threshold by the amount 

E x = ll/{2^). (15) 

We can therefore interpret this expression as the binding energy of the resonance [19]. 

The observed decays of the X imply that there are inelastic scattering channels, so a has 
a negative imaginary part. It can be parameterized in terms of the real and imaginary parts 
of 1/a as in Eq. (13). Our interpretation of X as a bound state requires 7 re > 0. The energy 
difference in Eq. (6) puts an upper bound on 7 rc : 

7re < 40 MeV (90% C.L.). (16) 

The upper bound on the width in Eq. (2) puts an upper bound on the product of 7 re and 

7im- 

7rc7im < (33 MeV) 2 (90% C.L.). (17) 

There is also a lower bound on the width of the X from its decays into D D -n- and D°D 0, y, 
which both proceed through the decay of a constituent D*. These decays involve interesting 
interference effects, but the decay rates have smooth limits as the binding energy is tuned 
to [11]. In this limit, the constructive interference increases the decay rate by a factor of 
2, so the partial widths of X into D°D°7i° and D°D°-f add up to 2T[D*°]. The width of 
D*° has not been measured, but it can be deduced from other information about the decays 
of D*° and D* + . Using the total width of the D* + , its branching fraction into D + n°, and 
isospin symmetry, we can deduce the partial width of D*° into D n to be 42 ± 10 keV. 
The total width of the D*° can then be obtained by dividing by its branching fraction into 
D°7i°: T[D*°] = 68 ± 16 keV. The sum of the partial widths of X into D D tt and D D°~/ 
is therefore 136 ± 32 keV. The resulting lower bound on the product of 7 re and 7; m is 

7rc7im > (7 MeV) 2 (90%C.L.). (18) 

By combining this with the upper bound on 7 re in Eq. (16), we can infer that 7 im > 1 MeV. 



III. HADRONIC STATES WITH NEARBY THRESHOLDS 



The state of the X can be written schematically as 

7f 



yl/2 

X) = =^=- (\D°D*°) + \D*°D )) + z h 2 \H), (19) 
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where Z is the probability for the X to be in the D D*°/D*°D state and Zh is the prob- 
ability for the X to be in another hadronic state H . The hadronic states H in (19) could 
include charmonium states, other charm meson pairs such as D ± D* T , states consisting of 
a charmonium and a light hadron such as J/tpp and J/ipcu, etc. An expansion of \X) in 
terms of hadronic states can be valid only if there is an ultraviolet cutoff on the energy 
difference with respect to the D°D*° threshold. The probabilities Zh depend on that cutoff. 
Universality implies that as a — > oo, Z approaches to 1 and Z H scales as 1/a [13]. 

The small binding energy of X compared to the natural energy scale m^/{2p) = 10 
MeV associated with pion exchange implies resonant S-wave interactions in the D°D*° and 
D*°D° systems. If there are other hadronic channels whose thresholds differ from the D°D*° 
threshold by less than 10 MeV, there would be resonant S-wave interactions in those channels 
as well. In this case, it would be necessary to treat all the resonating channels as a coupled- 
channel system, with a large elastic scattering length for each channel and a large transition 
scattering length for each pair of channels. The hadronic states D ± D* T , J/ip p, and J/ipu 
have thresholds that are relatively close to the D°D*° threshold. The energy gaps between 



these other thresholds and the D°D*° threshold are 

m D ± + m D *T — (m D o + m D *o) = +8.1 ± 0.1 MeV, (20a) 

mj/^ + m p - (m D o + m D *o) = +1.4 ± 1.1 MeV, (20b) 

m w + m u) - (m D o + m D *o) = +8.2 ± 1.0 MeV. (20c) 



The small uncertainty in Eq. (20a) comes from using mass differences between charm mesons 
to calculate the energy gap. The uncertainties in Eqs. (20b) and (20c) are dominated by the 
uncertainty in 2mfl0. The energy gaps in the D ± D* =F and J/ipu channels are comparable to 
the natural energy scale of about 10 MeV associated with pion exchange. The energy gap in 
Eq. (20b) for the J/ip p channel is much smaller. However whether any of these channels can 
have resonant interactions with D°D*° or D*°D° is determined not only by the real parts of 
the energy gaps, which are given in Eqs. (20), but also by the imaginary parts, which can 
be obtained by replacing each mass m by m — iT/2, where T is the width of the particle. If 
there are large differences between the widths of the various particles, it is necessary only to 
take into account the largest width among the particles in each channel. The largest width 
in each of the three channels is 

rfD**] = 0.096 ± 0.022 MeV, (21a) 
T[p] = 150.3 + 1.6 MeV, (21b) 
I>] = 8.49 ± 0.08 MeV. (21c) 

For the D ± D* T and J/ipu channels, the magnitude |A| of the complex energy gap is com- 
parable to the natural energy scale 10 MeV associated with pion exchange between D and 
D*. The large width of the p makes | A| for the J ftp p channel much larger than the natural 
energy scale. 

Because the complex energy gap A for the other hadronic channels with nearby thresholds 
are comparable to or larger than the natural energy scale, these channels need not be taken 
into account explicitly in calculations of quantities that have nontrivial universal limits as 
a — > ±oo. Their dominant effects enter through the complex-valued scattering length a. A 
coupled-channel model that includes other hadronic states with nearby thresholds could still 
be useful for estimating nonuniversal quantities or for calculating nonuniversal corrections 
to the universal predictions. 
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FIG. 1: The probabilities and Z^ p for the J/ipu and J/ipp components of the X(3872) as a 
function of its binding energy Ex- The dots are the results from Swanson's model [14]. The dotted 

1/2 

curves have the scaling behavior E-£ and pass through the model result for Ex = 0.7 MeV. 



Swanson has constructed a model of the X and the hadronic states with nearby thresholds 
and used it to predict some of the properties of the X [14, 16]. In particular, he predicted 
correctly that the branching fraction for X — > <//■?/> 7r + 7T~7r° is comparable to that for X — > 
J /4) 7T+7T-. In addition to the channel D°D*° + D*°D°, Swanson's model includes D + D*~ + 
D* + D~, J/vpp, and J/ipu;. It includes the S-wave and D-wave channels for DD*, but 
only the S-wave channel for J/ipV, where V is the vector meson p or u. Thus the model 
has 6 coupled channels. The interactions between the hadrons are modeled by potentials: 
one-pion-exchange potentials for the S-wave and D-wave DD* channels and for transitions 
between those channels and Gaussian potentials for the transitions between the S-wave 
DD* channels and the S-wave J/ipV channels to simulate the effects of quark exchange. 
The one-pion-exchange potential is singular at short distances and it was regularized by an 
ultraviolet momentum cutoff A. The nonrelativistic Schrodinger equation for the 6 coupled 
channels was solved numerically. A bound state with the quantum numbers J PC = 1 ++ 
of X appeared when the ultraviolet cutoff exceeded the critical value A c = 1.45 GeV. The 
binding energy of the X could be adjusted by varying the ultraviolet cutoff. 

Swanson solved the coupled channel problem under the assumption that the p and u are 
stable hadrons with equal masses m p = m w = 782.6 MeV. The reason for using an unphysical 
value for m p is that the central PDG value from 2002 and earlier, m p = 771.1 MeV, is below 
the D°D*° threshold. If such a value had been used, it would have been necessary to treat 
J/ip p states as scattering states. This complication was avoided by using a value of m p 
above the D°D*° threshold. Note that in Eq. (20b), we have taken the updated 2004 PDG 
value m p = 775.8 ± 0.5 MeV [45], which gives a J/ipp threshold that is a few MeV above 
the mass of the X. 

Swanson calculated the probabilities for each component of the wavefunction of X for 
values of the ultraviolet cutoff that correspond to varying the binding energy E x from 0.7 
MeV to 23.2 MeV [14]. His results for the probabilities and Z^p are shown as dots in 
Fig. I. Since the binding energy of the X is known to be less than 1 MeV, only the lowest 
two values of Ex could be physically relevant. For the lowest value Ex = 0.7 MeV, the 
probabilities were Z^ = 9.6%, Z d ± d *t = 7.9%, and Z^, p = 0.86%. The total probability 
for the D°D*° and D*°D° components of the wavefunction is 81.6%. 

1 /2 

In Fig. 1, the dotted lines have the scaling behavior E^ predicted by universality and 
are normalized so that they pass through the dot at E x = 0.7 MeV. The probability Z^ 
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clearly exhibits the universal behavior. The probability Z^, p does not. This may be related 
to the fact that Z^ p is more than an order of magnitude smaller than Z^. Because of the 
weaker coupling of the X to isospin-1 states, the scaling region for isospin-1 states may not 
set in until a much smaller value of E X - 

Swanson estimated the partial widths for the decays of X into J / ip h, where h is the light 
hadronic state 7r + 7r~, 7r + 7r~7r°, 7r°7, or 7r + 7r~7, using a simple ad hoc recipe. The partial 
width into J/xjj h was taken to be the sum over the vector mesons V = p, oj of the product 
of the probability Z^y for the J/ipV component of the wavefunction and the partial width 
r[V — > h] for the decay of the vector meson: 



For the smallest value of the binding energy that was considered, E x = 0.7 MeV, the 
resulting estimates of the partial widths for decay into J/iph were 1290 keV, 720 keV, 
70 keV, and 13 keV for h = it + it~, 7r + 7r~7r°, 7r°7, and 7r + 7r~7, respectively. The ratio of 
the partial widths into J/ip 7r + 7r~7r° and J/tpn+n^ was predicted to be 0.56 for Ex = 0.7 
MeV. Remarkably, this prediction agrees with the subsequent measurement by the Belle 
collaboration given in Eq. (4) to within the experimental errors [35]. The approximately 
equal branching fractions are a fortuitous result of an amplitude for X — > J/ipu that is 
much larger than the amplitude for X — > J/tp p and an amplitude for oj — > 7r + 7r~7r° that is 
much smaller than that for p — > 7r + 7r~. The suppression of the amplitude for X — > J/ip p is 
related to the fact that in the isospin symmetric limit in which the mass difference between 
neutral and charged D's is neglected, there is binding in the isospin-0 channel but not in 
the isospin-1 channel [43]. 

Swanson has also used his model to calculate the rates for several other decay modes 
of the X [16]. The decay rate into J/ip'j has contributions from transitions to J/ipp and 
J/ipuj that can be calculated using vector meson dominance. It also has contributions from 
the annihilation of the u and u from the charm mesons that are the constituents of the X. 
Swanson's prediction for the partial width into J/ip'J for an X with a binding energy of 1 
MeV is 8 keV. Decay modes that receive contributions only from uu annihilation, such as 
i/)(2S) 7, KK*, and np, have much smaller partial widths. 

IV. DECAYS OF X INTO J/tp h 

In this section, we calculate the differential decay rates of the X into J/iph, where the 
hadronic system h is 7r + 7r~7r°, 7r + 7r~, 7T°7, or 7. We assume that these decays proceed 
through transitions of X to J/ip p and J/ipuj. We calculate the differential decay rates 
in terms of two unknown complex coupling constants using an effective lagrangian that 
reproduces the decays of the light vector mesons. 

A. Vector Meson Decay Amplitudes 

We assume that the decay of X into J/ip h, where h is a system of light hadrons, proceeds 
through transitions to J/ipV, where V is one of the vector mesons p or uj, followed by the 
decay of the vector meson into h. Because the mass of the X is so close to the threshold for 
J/ip V, the vector meson is almost on its mass shell. Any model that reproduces the decays 




(22) 



v 
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of the vector mesons should also accurately describe the decay of the virtual vector meson 
in the J/ipV component of X. In Ref. [46], the semileptonic branching fractions for the r 
lepton were calculated using an effective lagrangian for light pseudoscalar and vector mesons 
with U(3)xU (3) chiral symmetry. All the parameters in the effective lagrangian, aside from 
the pion decay constant, were determined directly from decays of the vector mesons p and 
uj. That same effective lagrangian can be used to calculate the partial widths of X into 
J /if) h. An updated determination of the parameters in that effective lagrangian is given in 
the Appendix. 

The T-matrix element for the decay of a vector meson V into the light hadronic state h 
can be expressed in the form 

T[V - h] = e»A,[V - h], (23) 

where ey is the polarization vector of the vector meson. The amplitude A^ for the decay 
p — > 7r + 7r~ is 

A^p -> 7T+7T-] = \ G V7T7T (p+ - . (24) 

The value of the coupling constant G V7rn is given in Eq. (Ala). The amplitude A^ for the 

decay uj — > 7r + 7r~7r° is 

+ - oi 4 v / 3(cos6» 1 , + v^sin^) v a p 

A„[oJ^n + n tt u ] = — e^ a(3 p u +P 1Po 

IT 

x (c v3n + G ^ C - F " (1 _ I [/p ( Sl2 ) + f p{s23 ) + f p ( S31 )})) , (25) 
where Si 2 , s 2 3, and s 31 are the invariant masses of the three different pion pairs and 

MS) = s-M v S +l m v T v (26) 

is a vector meson resonance factor that vanishes at s = 0. We have denoted the 4-momenta 
of 7r + , 7T~, and 7T° by p + , p_, and po, respectively. The pion decay constant is F n = 93 
MeV, the values of the parameters C v3w and G V7T7T C vvn F%/ml are given by Eqs. (Alb) and 
(Ale), and the value of the light vector meson mixing angle 9 V is given by Eq. (A3). The 
amplitudes A^ for the radiative decays of the vector mesons are 

AM - vr+7] = il (c my + G ^ C - F2 A e^Q^ (27a) 

or v \ Tfl v J 

. r n i 4(cos^ + ygsingje / g^C^g \ « 

7T 7J = ^=— I G vnl H I e^apQ p e^, (27b) 

where Q and p are the 4-momenta of the vector meson and the pion and e 7 is the polarization 
vector of the photon. The values of the parameters C W7 and G vl C vvn F%/ml are given by 
Eqs. (A2b) and (A2c). 

The amplitudes A v in Eqs. (24), (25), and (27) all satisfy Q V A U = 0, where Q is the 
4-momentum of the vector meson. This condition is satisfied in any model consistent with 
vector meson dominance. The assumption of vector meson dominance is that the amplitude 
for the production of a real photon in a hadronic process can be expressed as the sum of 
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over vector mesons V of the amplitude for producing V multiplied by a coupling constant 
for the transition V — > 7. The condition Q h ' A v = is required for the gauge invariance of 
the resulting amplitude for real photon production. 

The T- matrix element for X to decay into J/ip and a light hadronic system h through a 
virtual vector meson resonance V can be expressed as 

T[X Jm = AJX - mV] Q2 _-fl imvT MV - hi (28) 

where Q is the total 4-momentum of the hadronic system h or, equivalently, of the virtual 
vector meson. We have used the condition Q U A V = to simplify the numerator of the vector 
meson propagator. The quantum numbers of the particles, together with Lorentz invariance, 
constrains the amplitude for X — > J/ip V to be the sum of two terms. One of them is 

A,[X - J /i, V] = G X1PV e^pQ^^ef, (29) 

where ex and e$ are the polarization 4- vectors of the X and the J / ip and Gxtpv is a dimen- 
sionless constant. The contraction of this amplitude with the polarization vector e v of the 
vector meson reduces in the rest frame of the vector meson to Gxt(>v m v e x • (e^ x e v)* ■ The 
other independent amplitude A^ has the Lorentz structure s llvol pP v e x e^ . In the rest frame 
of the X, its contraction with e v is m X £x • ( e v x e v)*- Since the mass of the X is so close to 
the threshold for J/ipV, the rest frames of the X and V are essentially identical. Thus the 
two independent Lorentz structures are essentially equivalent for decays that are dominated 
by the vector meson resonance. They give similar predictions for the partial widths for X 
into J/iph for h = 7r + 7r~7r°, 7r + 7r~, or 7r°7. The amplitude in Eq. (29) has the advantage 
that it is also consistent with the constraint Q^A^ = required by vector meson dominance. 
Thus this amplitude can be used to calculate the decay of X into J/tp'y. We therefore take 
the transition amplitude for X into J/t/iV to be the expression in Eq. (29). 

B. Decay into J/i/>7r+7r - 

We assume that the decay of X into <//■?/> 7r + 7r~ proceeds through a transition of X to 
J/ijj p. The T-matrix element is then given in terms of the unknown coupling constant Gxtpp 
by Eqs. (28) and (29) with V = p. The expression for the amplitude A v for p — > 7r + 7r~ is 
given in Eq. (24). We obtain the decay rate by squaring the amplitude, summing over spins, 
and integrating over phase space. The differential decay rate into <//■?/> 7r + 7r~ as a function 
of the invariant mass Q of the two pions is 

dT lx j,, + \Gx^ P \ 2 Gl n (Q"-Aml)^X^(mx,m^Q) 

dQ [ ~^ 9216vr 3 m^m2 (Q2 - m 2) 2 + m 2 p p 

x [(m 2 x + m%)(m\ - m%f - 2(m x - ±m\m% + m%)Q 2 + {m\ + mJ)Q 4 ] ,(30) 
where X(x, y,z) is the triangle function: 

X(x, y, z)=x 4 + y 4 + z 4 - 2(x 2 y 2 + y 2 z 2 + z 2 x 2 ). (31) 
After integrating over the pion invariant mass, the decay rate is 

T[X -> J/ipn+n-} = |G x ^| 2 (223keV). (32) 
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Q(n n ) [GeV] 

FIG. 2: Invariant mass distributions for the pions in the decay X — > J/ipir + 7r~. 

The shape of the pion invariant mass distribution for the decay of X into J/ip7T + 7T~ 
is shown in Fig. 2. Its qualitative features are dominated by the phase space factor 
\ l l 2 {m Xl m^, Q), which cuts the distribution off at the endpoint Q = mx — m^, and the 
vector meson resonance factor, which has its maximum at Q — m p just outside the kinematic 
region. Most of the support for dT/dQ comes from within T p of the upper endpoint. 



C. Decay into J/ij) 7r+7r 7r° 

We assume that the decay of X into J/ip 7r + 7r~7r° proceeds through a transition of X 
to J/ipuj. The T-matrix element is then given in terms of the unknown coupling constant 
Gxipui by Eqs. (28) and (29) with V — uj. The expression for the amplitude for u — > 3n is 
given in Eq. (25). We obtain the decay rate by squaring the amplitude, summing over spins, 
and integrating over phase space. The differential decay rate into J/ifj7r + ir~TT° as a function 
of the invariant mass Q of the 3 pions can be reduced to a 2-dimensional integral: 



dQ [X 



IGx^jcosdy + y^sin^) 2 A 1/2 (m x , m^, Q) 
30727r 5 m^m 2 F% 



^ J 3072^m x mlFt Q[(Q* - m 2 ) 2 + m 2 r 2 ] 

x [(rn 2 x + ™l)(m 2 x - ml) 2 - 2(m x - ±m\m% + m%)Q 2 + {m\ + mj)g 4 ] 

x J ds 12 J ds 23 [su S23S31 - rn 2 n (Q 2 - m 2 ) 2 ] 



x 



CUr + (1 - I [ /p ( Sl2 ) + f p ( S23 ) + / p ( S31 )]) 



(33) 



where S12, S23, and S31 are the squares of the invariant masses of the three pairs of pions. We 
have suppressed the limits of integration in the integrals over s\ 2 and s 23 . After integrating 
over the pion invariant masses, the decay rate is 



r[X^ J/V'tt+tt-tt ] = |G x ^| 2 (19.4keV). 



(34) 



The shape of the pion invariant mass distributions for the decay of X into J/-^7r + 7r n° 
is shown in Fig. 3. Its qualitative features are dominated by the phase space factor 



11 



0.5 



0.6 




Q(7t n Jl ) [GeV] 

FIG. 3: Invariant mass distributions for the pions in the decays X — > J/ip ir + n~n . 

\ l / 2 {mx-,rn^^Q)^ which cuts the distribution off at the endpoint Q = m x — m^, and the 
vector meson resonance factor, which has its maximum at Q — just outside the kine- 
matic region. Most of the support for dT/dQ comes from within a few widths of the 
upper endpoint. 

The ratio of the decay rates in Eqs. (32) and (34) is 



T[X -> J/jj TT+TT-TT 01 

T[X -> J/4)n+n- 



0.0870 



\G 



Xi/>w\ 



\G 



(35) 



By comparing this to Belle's result in Eq. (4) for the ratio of the branching fractions, we 
can obtain an estimate of the ratio of the coupling constants: 



\G 



Xipw\ 



G 



« 11.5 ±5.7. 



Xi,p\ 



(36) 



D. Decay into «//?/> 7r°7 

We assume that the decay of X into J/ipn ^ proceeds through transitions of X to J/ip p 
and J/ipu. The T-matrix element is then given by Eq. (28) summed over V = p, uu. The 
amplitudes for V — > 7r°7 are given in Eqs. (27). The differential decay rate with respect to 
the invariant mass Q of the 7r°7 is 



dr «em (C CT7 + G vl C vv ^F 2 /m 2 v f (Q 2 - m 2 ) 3 A 1/2 (mx, ny, Q) 

dQ [ ^ /V,7F 7J ~ 648vr 2 m^m 2 F 2 Q 



x [(rr& + mj)(m^ - m 2 ) 2 - 2(m^ - 4m^m 2 + m%)Q 2 + (m^ + m 2 )Q 4 ] 



x 



G 



Q 2 - m 2 + im p r / 



+ 



Gxv^V^cos^ + v^sin^. 



Q 2 - m 2 + imj: u 



(37) 



After integrating over the 7r°7 invariant mass, the decay rate is 

r[x^j/vA] = [|G^I 2 + °- 026 \Gxipp\ 2 

+ (0.163 cos + 0.215 sin <p)\G x ^\\G x ^p\\ (3.24 keV), (38) 
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where exp(i<f)) is the relative phase between Gx-ipw and Gx-ipp- The estimate of the ratio 
\Gxipu)\ 2 /\Gxipp\ 2 in Eq. (36) suggests that the |Gx^| 2 term in Eq. (38) dominates. If this 
is the case, the branching fraction for the decay of X into J/-07r°7 should be smaller than 
that for J/tl)ir + ir~iT by a factor of about 0.17. 



E. Decay into J/tjij 



Having chosen the transition amplitude in Eq. (29) so that it satisfies Q^A^ = 0, we can 
use vector meson dominance to calculate the partial width for the decay of X into J/ipl- 
The T-matrix element is 



G F 2 e 



G 



rat 



im p r p 



+ 



G 



cos 



m. 



L w 



(39) 



where Q is the 4-momentum of the photon. The value of the coupling constant G VJ is given 
in Eq. (A2a). The result for the decay rate is 



2 \3 



a em G 2 VJ F*(m x + m%){m x - ra\) 



24ra x ra 2 b 



x 



G 



X^p 



rat 



ira p T p 



+ 



G 



Xipui 



COS 



e v /y/3 



rat 



If the widths in the vector meson propagators are neglected and if we use m p 
decay rate in Eq. (40) reduces to 

V[X -> J/^ 7 ] = |Gx^ + 0.30Gx^| 2 (5.51keV). 



(40) 
m U) the 

(41) 



Our estimate in Eq. (36) implies that is much larger than \Gx%j} P \- However, the 

larger magnitude of Gxi/ju is compensated by the vector meson mixing factor cos9 v /\/3 = 
0.30, so the Gxt(,p and Gx^ui terms may be equally important. Using the partial widths in 
Eqs. (32) and (34), we can relate the branching fractions for J/-07 to those for J/-07T + 7r~ 
and J/ip 7r + 7r~7r°: 

Br[X -> J/^ 7 ] = 0.025 Br[X -> J/ijj tt+tt-] + 0.026 Br[X -> J/^ tt+tt-tt ] 

+0.050 cos (Br[X -> J/^ tt+tt-] Br [X -> J/0 tt+tt-tt ]) 1/2 , (42) 

where exp(i0) is the relative phase between Gx^d and Gx-ipp- This prediction is compatible 
with the measurements of the branching ratios in Eqs. (4) and (5) if the angle is small. 



F. Factorization of short-distance decay rates 

The decay modes of the X(3872) can be classified into long-distance decays and short- 
distance decays. The long-distance decay modes are 7J°Z) 7r and D°D 0/ -f, which proceed 
through the decay of a constituent D*° or D*°. These decays are dominated by a component 
of the wavefunction of the X in which the separation of the D and D* is of order \a\. 
These long-distance decays involve interesting interference effects between the D°D*° and 
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D*°D° components of the wavefunction [11]. The short-distance decays are dominated by a 
component of the wavefunction in which the separation of the D and D* is of order or 
smaller. Examples are the observed decay modes J/tp 7i + n~, J/ip 7r + 7r _ 7r°, and J/ipj. 

Short-distance decays of the X into a hadronic final state H involve well-separated mo- 
mentum scales. The DD* wavefunction of the X involves the momentum scale l/|a| set by 
the large scattering length. The transition of the DD* to H involves momentum scales 
and larger. The separation of scales |a| l/m n can be exploited by using a factorization 
formula for the decay rate [19]. In limit \a\ ^> l/m n , the leading term in the T-matrix ele- 
ment for the decay X — > H can be separated into a short-distance factor and a long-distance 
factor: 

T[X -> i?] = A short [X ^H]x A x . (43) 

The short-distance factor -4 s hort in Eq. (43) has a well-behaved limit as \a\ — > oo. The long- 
distance factor Ax is the universal amplitude given in Eq. (12). If the complex scattering 
length is parameterized as in Eq. (13), this factor is 

A X = (fht/v) (7re + *7im) 1/2 ■ (44) 

When applied to decays of X into J/ip and light hadrons, the factorization formula in 
Eq. (43) implies that the coupling constants Gx^ P and Gx^w have a long-distance factor 
Ax- 

The factorization formula for the T-matrix element in Eq. (43) implies a factorization 
formula for the decay rate: 

T[X -+H] = r short [X -> H\ x \A X \ 2 . (45) 

The short-distance factor r s h or t in Eq. (45) has a well-behaved limit as |a| — > oo. Using the 
expressions in Eqs. (15) and (14) for the binding energy and the width of the molecule, the 
long-distance factor in Eq. (45) can be expressed as 

\A X \ 2 = [Ex + T x /(16Ex)} 1/2 . (46) 

Predictions for the rates for short-distance decays of the X can be obtained from models 
for low-energy hadrons in which the parameters have been tuned to obtain a small bind- 
ing energy Ex, such as Swanson's model [14]. In such models, calculations using the most 
straightforward numerical methods tend to become increasingly unstable as the binding en- 
ergy is tuned toward 0, because the small binding energy results from a delicate cancellation. 
The factorization formula in Eqs. (45) and (46) can be useful for extrapolating the predic- 
tions of a model to other values of the binding energy Ex- In many models, it is difficult to 
take into account effects of the width Fx of the molecule. Given the prediction of a model 
in which the width has been neglected, the factorization formula in Eqs. (45) and (46) can 
be used to take into account the nonzero width Tx consistently. 

In order to use the factorization formula in Eqs. (45) and (46) to extrapolate a partial 
width calculated using a model to other values of the binding energy and the width, the 
calculation must be carried out for small enough binding energy that the model is in the uni- 
versal scaling regime where observables scale as powers of the binding energy. For example, 
the probabilities for components of the wavefunction other than D°D*° and D*°D° should 
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scale as E-^ . In Swanson's model with Ex = 0.7 MeV, this universal scaling behavior is 
satisfied by the probability Z^ w but not by Z^ p , as is evident in Fig. 1. The delayed onset 
of the universal behavior for the probability perhaps be attributed to the weaker 

coupling of X to isospin-1 states. In the next section, we will use Swanson's result for 
to estimate the coupling constant Gx^w 

V. PARTIAL WIDTH FOR X -» J /if; tv+tv'tv 

The partial widths of the X calculated in Section IV are expressed in terms of unknown 
coupling constants Gx-ipp and Gx^w In this section, we use a simple 2-channel scattering 
model to show that | Gx^w | can be deduced from the probability Z^ for the J/ipu component 
of the X. We then use the probability Z,^ in Swanson's model to give a quantitative 
prediction for the partial width for the X to decay into J/ipn + n^n . 



A. Two-channel scattering model 



Cohen, Gelman, and van Kolck have constructed a renormalizable effective field theory 
that describes two scattering channels with S-wave contact interactions [47]. We will refer 
to this model as the two- channel scattering model. An essentially equivalent model has been 
used to describe the effects of A A states on the two-nucleon system [48] . The parameters of 
this model can be tuned to produce a large scattering length in the lower energy channel. It 
can be used as a simple model for the effects on the D D*°/D*°D system of other hadronic 
channels with nearby thresholds, such as J/ip p and J/ipoj. 

The two-channel model of Ref. [47] describes two scattering channels with S-wave contact 
interactions only. We label the particles in the first channel la and 16 and those in the 
second channel 2a and 26. We denote the reduced masses in the two channels by /x and \x 2 . 
Renormalized observables in the 2-body sector are expressed in terms of 4 parameters: three 
interaction parameters an, a22, and ai2 = 021 with dimensions of length and the energy gap 
A between the two scattering channels, which is determined by the masses of the particles: 



A = m 2a + m 2 b - (m la + m lh ). 



(47) 



The scattering parameters in Ref. [47] were defined in such a way that a n and a 2 2 reduce in 
the limit a 12 — > ±00 to the scattering lengths for the two channels. The truncated connected 
transition amplitude A(E) for this coupled-channel system is a 2 x 2 matrix that depends 
on the energy E in the center-of-mass frame. If that energy is measured relative to the 
threshold m\ a + mu, for the first scattering channel, the inverse of the matrix A(E) is 

1/on + y/-2fJ,E] y//i^2~/ai2 

y/mpz/au A* 2 [ — l/a 2 2 + A/2/i 2 (A - E) ] 

The square roots are defined for negative real arguments by the prescription E ■ 
with e — > + . The explicit expressions for the 11 and 12 entries of this matrix are 

l/a 22 + v /2 f i 2 {A- E) 1 



(48) 



re 



* 2vr / 1 , 1 

A U (E) = — + y/-2(j,E - — 2~ 



A 12 (E) 



2ir 



a n 
1 

Ol2 



12 



Ol2 



an 



+ y/-2flE 



a 2 2 



+ y/2fl 2 (A - E) 



(49a) 
(49b) 
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The amplitudes denned by (48) are for transitions between states with the standard non- 
relativistic normalizations. The transitions between states with the standard relativistic 
normalizations are obtained by multiplying by a factor ^/2mj for every particle in the ini- 
tial and final state. The T-matrix element T n (p) for the elastic scattering of particles in 
the first channel with relative momentum p is obtained by evaluating An(E) at the en- 
ergy E = p 2 /(2fi). The scattering length is determined by the T-matrix element at p = 0: 
^ii(O) = —2na/fi. The inverse scattering length 1/a is therefore 



1 1 1 

a an O'u 



\PM^K- l/a 22 l . (50) 



If the matrix A(E) given by Eq. (48) has a pole on the physical sheet at E = — k 2 /(2/x), 
there is a bound state below the scattering threshold for the first channel with binding 
energy Ex = k 2 /(2ji). The binding momentum k satisfies 



k = h — 

an a 12 



- l/a 22 + + Ml*)* 2 ] • (51) 



The momentum-space wavefunction ip{p) for the bound state is a column vector whose 
two components are the amplitudes for the bound state to consist of particles with relative 
momentum p in the first and second channel, respectively. The wavefunction can be deduced 
from the behavior of A(E) near the bound-state pole: 



(52) 



The components Axi and .4x2 of the column vector are the amplitudes for transitions from 
the bound state to particles in the first and second channels, respectively. They satisfy 

(x[-l/a u + k] Axi + WW21 012] -4x2 = 0. (53) 

Because the only interactions in the two-channel model are contact interactions, the depen- 
dence of the wavefunction on the relative momentum of the constituents comes only from 
propagators. The wavefunction can be expressed in the form 

- n {2^A x ^:72i^]-) • (54) 

where iV is a normalization constant. The normalization condition 

/^(I^(P)| 2 + I^2(P)| 2 ) = 1 (55) 

can be expressed as Z\ + Z2 — 1, where Zi and Z2 are the probabilities for the bound state 
to consist of the particles in the first and second channels, respectively. The probability Z x 
for the first channel is given by 

1 V^+(^ 
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B. Two-channel model with large scattering length 



In the two-channel model of Ref. [47], a large scattering length a in the first channel can 
be obtained by fine-tuning the parameters an, a 2 2, 012, and A. The natural momentum scale 
A associated with low-energy elastic scattering in the first channel is set by the magnitudes 
of a,ii , a i2i an d (2/U2A) 1 / 2 . There are various ways to tune the parameters so that \a\ 
is large compared to A -1 . For example, a can be tuned to ±oo by tuning the scattering 
parameter an to the critical value — a^tv^A^A — 1/022]- 

As a is tuned to be much larger than the natural momentum scale, the amplitude Au(E) 
for \E\ <C A 2 /(2/i) approaches the universal expression given in Eq. (9). The solution to 
Eq. (51) for the binding momentum k approaches 1/a, so if a > 0, there is a bound state 
with the universal binding energy in Eq. (10). The first component of the wavefunction in 
Eq. (54) approaches the universal expression in Eq. (11), while the probability of the second 
component approaches as 1/a. The amplitude Axi for the transition from the bound state 
to particles in the first channel also approaches the universal amplitude Ax in Eq. (12). 

There are also universal features associated with transitions from the bound state to 
particles in the second channel. If \a\ 3> A -1 , the leading term in the amplitude for the 
transition of the weakly-bound state X to particles in the second channel is 

A X2 = - [ v /^a _ i/ a22 1 ' 1 Ax} (57) 



Ol2 



where Ax is the universal amplitude given in Eq. (12). This equation is a factorization 
formula that expresses the transition amplitude as the product of a short-distance factor 
and the universal long-distance factor Ax- Using Eq. (56), the probability Z 2 = 1 — Z 1 for 
the bound state to consist of particles in the second channel reduces to 



a 



W/i) 

2 2V / 27^A~ 



v / 2a^A- l/a 22 l 2 -- (58) 



Note that the probability Z 2 differs from |«4x2| om y by kinematic factors: 



\A 



X2\ 



|2 



y / 8^ 2 A/^Z 2 . (59) 



This relation also follows directly from the wavefunction in Eq. (54) if we use the fact that the 
normalization factor N approaches 1 as a — > 00. Thus the relation between the probability 
and the transition amplitude in Eq. (59) is not specific to the 2-channel model. It applies 
more generally to any 2-particle component of the bound state whose wavefunction can be 
approximated by {p 2 + 2// 2 A) _1 , where A is the energy gap. It requires only that A is small 
enough that the interaction in that channel can be approximated by an S-wave contact 
interaction at momenta comparable to v^/^A. 



C. Partial width into J/tj) 7v + tv 7r° 

We can use results from Swanson's model to estimate \Gx^uj\, thereby determining the 
unknown constant in the expression in Eq. (34) for the partial width for X — > J/tp 7r + 7r~7r°. 
The relativistic amplitude for the transition from X to J/ipu is given by the contraction 
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of the amplitude A^ in Eq. (29) with a polarization vector for the uj. The corresponding 
nonrelativistic amplitude Ax-ipw is the analog of the transition amplitude Ax2 in Eq. (57) 
for the 2-channel model. In the rest frame of the X, the relativistic amplitude differs from 
the nonrelativistic amplitude by a factor of \/2vrii for every external particle: 

e^A^X -> J/^uj] = {%m x m^m w ) xl2 Ax^- (60) 

Using the expression for the amplitude A p in Eq. (29) and the fact that the rest frame of 
the X is almost identical to that of the uj, the left side of Eq. (60) is 

tfA^X -> J/ipoj] = G x ^m u e x • (ey, x e w )*. (61) 

The transition amplitude Ax^ on the right side of Eq. (60) must have the same dependence 
on the polarization vectors of the J/ip and uj. There are two independent pairs of spin states 
for J/ip and uj that couple to any given spin state of X. If the analog of the factorization 
formula in Eq. (59) is summed over the spin states of the J ftp and uj, it gives 



J2 \Axi*,\ 2 = y^A^/z^V,, (62) 



spins 



where is the energy gap in Eq. (20c), ji^ is the reduced mass of the J/ip and uj, and 
is the probability for the J/ipu component of X. Squaring both sides of Eq. (60) and 
summing over the spin states of J/ip and uj, we get 



2mJ\G x ^\ = Vo-Km x {m^ + m aJ )y'2A^//i^ Z^. (63) 

Inserting Swanson's result Z^ = 9.6% for Ex = 0.7 MeV and using the factorization 
formula in Eqs. (45) and (46), this reduces to 



Ex + T x /(16E X )\ 



1/2 



10^ = 9.59^ Q - y - j . (64) 

Inserting the result into the expression in Eq. (34), we get a quantitative result for the partial 
width: 

r[X J hl , = (222 keV) ( ^M ) " 2 . (65) 

We can use the result in Eq. (65) to set a lower bound on the partial width into 
J/ijj7r + 7T~7T . As a function of the binding energy Ex, the right side of Eq. (65) is min- 
imized at E x = T x /4. The lower bound on the width is T x > 2T[D*°] = 136 ± 32 keV. 
Thus the lower bound on the partial width into J/ip 7r + 7r~7r° in Swanson's model is about 
58 keV. 

As is evident in Fig. 1, Swanson did not calculate the probability Z^p for the J/ipp 
component of X for a binding energy small enough to be in the scaling region where Z^ p 

1 /2 

scales like E x . If he had, we could use an equation analogous to Eq. (63) to determine 
\Gxipp\- If we assume that the smallest binding energy considered by Swanson is close to 
the scaling region, we can use his value Z^ p = 0.86% for E x = 0.7 MeV to estimate \Gx^ P \- 
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In the analog of Eq. (63), we should set A^ p = rather than using the value A^ p in 
Eq. (20b), because Swanson set m p = in his calculation. The resulting estimate is 

\Gx+„\ -0.86^ J . (66) 

We can insert this estimate into Eq. (32) to get an estimate of the partial width for decay 
into J ftp 7r + 7r~. We can also insert this estimate of \G x ^p\ 2 and the value of \G X ipu,\ 2 from 
Eq. (64) into Eqs. (38) and (41) to get ranges of estimates of the partial widths for the 
decays into <//■?/> 7r°7 and J/V'T- The ranges arise from the unknown relative phase between 
Gxipui and Gxipp- 



VI. SUMMARY 

Evidence is accumulating that the X(3872) is a loosely-bound S-wave molecule corre- 
sponding to a C = + superposition of D°D*° and D*°D°. Because its binding energy is 
small compared to the natural energy scale associated with pion exchange, this molecule has 
universal properties that are completely determined by the large scattering length a in the 
C = + channel of D°D*° and D*°D°. 

We have analyzed the decays of X into J/ip plus light hadrons under the assumption 
that X is a DD* molecule and that these decays proceed through transitions to J/ip p and 
J/i/juj. The differential decay rates were calculated in terms of unknown coupling constants 
Gxt(,p and Gx^ by using an effective lagrangian that reproduces the decays of the light 
vector mesons. The dependence on the unknown coupling constants enters only through 
multiplicative factors, so the angular distributions are completely determined. 

Quantitative predictions of the partial widths for the decays of X into J/ip plus light 
hadrons require numerical values for the coupling constants Gx^p and G^. We pointed 
out that the dependence of these coupling constants on the binding energy E x and the 
total width T x are determined by factorization formulas. We showed how \Gx^uj\ 2 could be 
determined from the probability for the J/ipu component of X in Swanson's model. We 
used this result to give a quantitative prediction for the partial width for X — > J/ip7i + 7i^7i 
as a function of Ex and Fx- 



APPENDIX A: DECAY AMPLITUDES FOR VECTOR MESONS 

In this appendix, we present an updated determination of the coupling constants in the 
effective lagrangian for the light pseudoscalar and vector mesons that was used in Ref. [46] to 
calculate the semileptonic branching fractions for the r lepton. The same effective lagrangian 
is used in Section IV to calculate the decay rates of the X into J /if) and light hadrons. 

The pion decay constant = 93 MeV and the hadron masses have all been determined 
accurately [45]. The other parameters in the effective lagrangian can be determined from the 
partial widths for decays of p°, p^, and u given in Table I. The most useful combinations of 
the parameters in the amplitudes for the decays of the vector mesons into pions in Eqs. (24) 
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Decay mode 


Partial width 




1 o i 1 ft MpV 


jJ ? O O 




F ' " 1 


67 6 + 7 5 keV 

\J 1 • U _1_ J . v7 JA.C V 


uj — ► e + e~ 


0.60 ± 0.02 keV 


w — > 7T°7 


0.76 ± 0.05 MeV 




0.82 ± 0.20 keV 


a; — > 7T + 7r _ 7r° 


7.56 ± 0.093 MeV 



TABLE I: Inputs that are used to determine the coupling constants in the vector meson decay 
amplitudes. The partial widths are taken from Ref. [45]. 



and (25) are 

G V7vn = 11.99 ±0.06, (Ala) 

C v37V + G VW7V C vvw F 2 /m 2 v = (8.03 ± 0.48)/(16vr 2 ), (Alb) 

G V7vn C VV7V F 2 /m 2 v = (10.2±1.3)/(16tt 2 ). (Ale) 

The coupling constant G V1 associated with vector meson dominance and the most useful 
combinations of parameters in the amplitudes for the radiative decays of the vector mesons 
in Eqs. (27) are 

G V1 = 14.01 ±0.11, (A2a) 
C V7Ty + G vy C vvw F 2 /m 2 v = (7.99±0.45)/(16vr 2 ), (A2b) 
G vl C VV7V F 2 /m 2 v = (11.9±1.5)/(16tt 2 ). (A2c) 

The vector meson mixing angle is given by 1 

cos 9 V = 0.51 ±0.01. (A3) 

Another function of 9 V that is often encountered is cos 9 V + V2 sm9 v ~ 1.73 ± 0.01. The 
errors in the parameters in Eqs. (Al), (A2), and (A3) are determined using the uncertainties 
in the measurements of the vector meson decay widths only. The uncertainties in the 
hadron masses and the pion decay constant are negligible in comparison. Variations in the 
parameters associated with U(3) x U (3) symmetry breaking are neglected in this analysis. 

The inputs that were used to determine the parameters in Eqs. (Al), (A2), and (A3) are 
listed in Table I. Following Ref. [46] , we determine the parameters by the following steps: 

1. The coupling constant G V7Tn in Eq. (Ala) is determined from the partial width for 

P — > 7T + 7T~: 

T[p - vr+vr-] = (1 - 4m 2 /m 2 ) 3/2 . (A4) 

2. The coupling constant G V1 in Eq. (A2a) is determined from the partial width for 

p — > e + e~: 

4vra 2 G 2 F 4 

T\p^e+e-]= -JJ \ (A5) 



1 The cosine of the angle 6 V here is the sine of the vector meson mixing angle used in Ref. [46] . 
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3. The combination of parameters in Eq. (A2b) is determined from the partial width for 

p~ — > 7r~7: 

T[p -> TT 7] = 27F / (^CW 7 + — ^ J (1 - . (A6) 

4. The combination of parameters G V1 C VV7T F 2 /ml in Eq. (A2c) is determined from the 
ratio of the partial widths for uj — > and uj — > 7r°7. The possibility of a 
relative phase between C vn7 and G V1 C VV7T F 2 /m1 is ignored. The partial width for 
uj — > 7r°/! + /!~ is 

r> - AV] = / ^ / (faasEI^ - AV1I 2 (A7) 

The squared amplitude, averaged over initial spin states and summed over final spin 
states, is 

J2\A[u -> ttWII 2 = 128 9 ^^ (cosfl„ + v^sin^) 2 

x [(s 2 3 + 4mJ) ((m 2 - s 23 - m 2 n ) 2 - 4to 2 s 23 ) + s 23 (s 12 - s 31 ) 2 ] 

(A8) 



1 

x 



s 2 



C C F 2 



c 7 + ° 7 7 

m 2 



where Si 2 , s 2 3, and s 31 are the squares of the invariant masses of the 7r°/i + , \i , and 
fi~ir°, respectively. The partial width for uj — > 7r°7 is 

I> - A] = 3(cos6>„ + ^sin^) ^^:^^ 3 3 r[p- - vr-7], (A9) 

where r[p~ — > 71-^7] is given in Eq. (A6). Note that the factor (cos0„ + v^sin^) 2 
cancels in the ratio of Eqs. (A7) and (A9). 

5. The combination of parameters G vn7V C vvn F 2 jm\ appearing in Eq. (Ale) is determined 
by multiplying the combination of parameters in Eq. (A2c) by the ratio G V7T7T /G V1 
obtained from Eqs. (Ala) and (A2a). 

6. The combination of parameters in Eq. (Alb) is determined from the ratio of the partial 
widths for uj to decay into 7r + 7r~7r° and 7r°7 and from the value of the combination 
of parameters in Eq. (Ale). The possibility of a relative phase between C„ 37r and 
G V7T7V C VV7T F 2 /ml is ignored. The partial width for uj — > 7r°7 is given in Eq. (A9). The 
partial width for uj — > 7r + 7r~7r° is 

r[uj - TtVtT-] = 256 ^ 3m 3 J ds 12 J dS2sZ\A[u - TT + TT-TT ]! 2 . (A10) 

The squared amplitude, averaged over the spin states of uj, is 

^,, r + - 0112 4(cosfl„ + y^sin^) 2 . 2 2 2 2 , 



x 



2 



+ (1 _ 1 [/p(si2 ) + /p ( S23 ) + /p(s3l)] ) 



m v 



2 



(All) 



Note that the factor (cos6*„ + v^sin^) 2 cancels in the ratio of Eqs. (A10) and (A9). 
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7. Finally, the cosine of the vector meson mixing angle in Eq. (A3) is determined from 
the ratio of the partial widths for cu — > e + e~ and p° — > e + e~: 

r> _> e + e-] = C ° s2 ^ m ' r[p° -> e + e-]. (A12) 
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